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In their comment on the paper (Phys. Rev. B 65, 153403 (2002); cond-mat/0110154), Alexandrov 
and Bratkovsky (cond-mat/0207173) argue that they correctly took into account the chemical poten- 
tial oscillations in their analytical theory of combination frequencies in multiband low-dimensional 
metals by expanding the free energy in powers of the chemical potential oscillations. In this reply, 
we show that this claim contradicts their original paper (Phys. Rev. B 63, 033105 (2001)). We 
demonstrate that the condition given for the expansion is mathematically incorrect. The correct 
condition allows to understand the limits of validity of the analytical theory. 



It is well-known that, due to the quantization of the electronic energy spectrum of metals into discrete Landau 
levels in the presence of a magnetic field, the chemical potential is expected to oscillate with the magnetic field when 
the number of electrons is kept constant. In two-dimensional (2D) multiband metals, it has been predicted initially in 
CO the framework of numerical studies 1 ' 2 that such chemical potential oscillations are responsible for the appearance of 
additional oscillations of the magnetization whose frequencies are combinations of the independent band frequencies 
fa- The number of electrons is fixed experimentally, independently of the dimensionality of the energy spectrum. For 
this reason, combination frequencies are in principle possible as well in 3D metals as in 2D multiband metals. However, 
chemical potential oscillations effects are not observed in multiband 3D metals. So, the mechanism of appearance of 
combination frequencies has to clearly point out the difference between low-dimensional metals and 3D metals. 
Alexandrov and Bratkovsky are the first authors to propose an analytical derivation for these combination fre- 
O , quencies 3 . However, they did not mention clearly the difference between low-dimensional metals and 3D metals in 
the framework of their theory. More precisely, an important point noticed for one-band 2D metals was missing in 
their derivation 3 : the analytical treatment of the chemical potential oscillations effects leads to a system of nonlinear 
equations 4,5 . In 3D metals, the resolutio n of these equations is trivial because the oscillating part jj, of the chemical 
potential is of the order of hio c ^J frw c j 1 ef (where co c is the cyclotron pulsation and Ef the Fermi energy). As a result, 
the magnetization oscillations are not sensitive to the chemical potential oscillations effects: the expression for the 
magnetization is (with high accuracy) the same for a fixed chemical potential p, as for a fixed number of electrons N. 
In 2D metals p ~ huj c , that leads to the observable difference between low-dimensional metals and 3D metals. This 
• fact has not been noticed in the article 3 and has motivated the analytical work 6 and its principal statement that the 
chemical potential oscillations appearing in the arguments of the Fourier components (of the grand canonical potential 
"j***; ■ or of the magnetization oscillations) were not taken into account by Alexandrov and Bratkovsky. 

In their comment 7 , these authors have addressed principally two criticisms to the paper 6 . First, (i) the chemical 
potential oscillations were correctly taken into account in their previous work 3 : they argue that they made in fact an 
expansion (as in the commented paper 6 ), but did not judge to mention explicitly this technical point 7 . Secondly, (ii) 
the analytical formula for the combination frequencies amplitude they derived is accurate even at zero temperature 
in a clean two-band 2D metal. We will show in this reply that the argument (i) is in contradiction with the original 
q \ paper 3 and is thus not receivable. The point (ii) is far from being correct. 
L" ' The main point of the Ref. 3 is to express the oscillating part F of the free energy as 
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where p is the total density of states of a two-dimensional multiband metal, and Q(p) is the oscillating part of the 
grand canonical potential, which is an explicit function of the chemical potential p. At a constant number of electrons 
N, the chemical potential p oscillates with the magnetic field and can be written as the sum of a constant part po 
(independent of the magnetic field) and an oscillating part p given by 
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As stressed in the Ref. 6 , this equation (2) is a self-consistent nonlinear equation to solve in order to determine the 
dependence of the chemical potential p on the magnetic field H. 

The grand canonical potential fl(p) depends on the chemical potential p through the expression 
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where 
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are the arguments entering in the Fourier components. Here uo ca is the cyclotron pulsation with the effective mass 
m a , A a is the a band-edge, and 

A r a = rn^ RT{r)RD{r) 



is the amplitude for the harmonic r in the band a with 
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T the temperature, \ r = 2ir 2 rkBT/fkd ca , and T the relaxation rate at H = 0. Since for a constant number of 
electrons the chemical potential oscillates with the magnetic field, so does the quantity f a . For this reason, the 
explicit expression (3) can not be seen as a Fourier series (f a is not a frequency, which by definition has to be 
independent of the magnetic field) . 

This difficulty can be overcomed by expanding the free energy F(p) (expression (1)) in powers of the oscillating 
part p of the chemical potential. By developing separately the first term and the last term of the right-hand side of 
the equation (1), we straightforwardly obtain keeping terms up to the order of pp 2 
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Therefore the free energy becomes 
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Obviously after equations (4) and (5), Fourier harmonics with combination harmonics are produced by both terms of 
the right-hand side of the expression (1), the contribution of the first term being reduced to one half due to a partial 
cancellation with the second term. 

Alexandrov and Bratkovsky argue 7 that the free energy F is expanded in powers of p -C Mo in their original work 3 ; 
in fact, the authors asserted just after giving the expression (1) (their Eq. (12)) in the Ref. 3 : "It is the last term, 
which yields combination Fourier harmonics with the combination frequencies f = rf a ± r' f a " ■ Thus, this sentence 
directly contradicts the fact that the free energy has been expanded in the paper 3 : the chemical potential oscillations 
were explicitly not taken into account in the trigonometric arguments. 

It is worth noting that a similar claim has been addressed independently by Kishigi and Hasegawa 8 , who wrote 
concerning the authors of Ref. 3 : "their result for the free energy [...] is formally correct but they did not take account 
of the magnetic- field dependence of the f a , which cannot be neglected in two-dimensional systems. As a result their 
analysis of the de Haas-van Alphen oscillation for the fixed N system (canonical ensemble) is insufficient and their 
conclusions on Fourier-transform intensities are incorrect" . 

The similitude of the formulae obtained in Ref. 6 and Ref. 3 has a simple explanation. We note by comparing the 
expressions (1) and (6) that the squared term changes sign after the expansion. It implies that if we ignore by hand 
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the oscillations of the chemical potential in the quantity f a , the amplitude obtained for the combination harmonics 
is fortuitously the same (a sign apart) after and before the development. 

The second point (ii) of discordance is in fact related to the validity of the analytical expansion of F (or of f2, or 
of the magnetization oscillations M) in powers of the oscillating part jl of the chemical potential /i. The condition of 
validity for it, \p,\ <C fio (which is fulfilled in the whole regime of magnetic quantum oscillations) given in 7 , is actually 
mathematically incorrect. The correct condition for performing the developments (4-6) is revealed when considering 
the explicit form (3) for CI(li). The basic point is that the oscillating part jl enters in the arguments of trigonometric 
functions for which the approximation 

cos(xo + x) w cos(aio) 

is valid under the condition \x\ <C 1 (and not \x\ -C xq, take e.g. 7r/2 <C x ). The expansion (6) is therefore valid 
provided that 

« 1 (7) 

for all significant harmonics r. Obviously, this latter condition (7) is much more restrictive than the condition given 
by the authors of the Ref. 7 and depends on the degree r of the harmonic (it is stronger for higher harmonics). It holds 
simultaneously for the expansion of all quantities which are functions of n/oJ ca through trigonometric arguments such 
as f2, the magnetization oscillations M or even jl itself (see equation (2)). 

It has already been noticed in Ref. 6 that jl is naturally reduced by the presence of multiple bands independently 
of the temperature or the impurity reduction factors. However, at very low temperatures and in clean 2D two-band 
metals (the most unfavorable case), many harmonics r are significant, and the condition (7) is not fulfilled for all terms 
because in this case |/2| is of the order of a few tenths of hw ca . Then, higher powers of jl have to be considered and 
the analytical treatment of the chemical potential oscillations effects is not obvious. In this regime of so-called strong 
chemical potential oscillations the validity of the analysis of the oscillations in terms of Fourier series is questioned at 
small but finite temperatures (this is not a property established a priori; the use of the Fourier analysis has thus to 
be justified). 

The importance of higher powers of jl has been demonstrated analytically and numerically in the case of one-band 
2D metals 4,5 . In these papers 4,5 the full nonlinear equation (2) is considered and solved at zero temperature. As a 
result, the drops of the magnetization oscillations occur at integer values of the ratio ^o/Hu c , while they occur at 
half-integer values when neglecting jl in the trigonometric arguments. Furthermore, harmonic amplitudes are found 
to differ strongly in the two situations, especially for high harmonics 5 . 

In the two-band 2D metals, the combination frequencies appear in the first order expansion in powers of jl. Alexan- 
drov and Bratkovsky argue 7 with the help of an analytical and a numerical estimates that higher powers of jl can be 
neglected to describe the magnetization oscillations with a relatively good accuracy even at zero temperature and in 
clean samples. In their analytical estimate, they only keep the first harmonic in the self-consistent equation (2) for 
the chemical potential oscillations, which is valid in fact at zero temperature only at a small Dingle factor i?£>(l). 
In this regime, the oscillations of jl are small and particularly smooth: for this reason the linear approximation is 
naturally expected to be quite good since the condition (7) is then fulfilled. On the contrary, in their numerical study 
they consider all the harmonics in Eq. (2), which is relevant in the regime of strong oscillations of jl. We can however 
cast some doubts on their result shown in the Fig. 1 of Ref. 7 : surprisingly, at small Dingle factors (i.e. when the 
expansion is possible), the accuracy of the linear approximation is increasing when the Dingle factor i?c(l) decreases; 
this is completely opposite to that can be expected. 

In conclusion, we have demonstrated that the chemical potential oscillations were not correctly taken into account 
in the Ref. 3 . Careful numerical studies are still needed to analyze the regime of strong chemical potential oscillations 
where the expansion of the quantities in powers of jl is not convergent. 

I thank V.P. Mineev for his advices and careful reading of the manuscript. 
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